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Abstract 

We study numerically the critical properties of the U(l)-Higgs lat- 
tice model, with fixed Higgs modulus, in the region of small gauge 
coupling where the Higgs and Confining phases merge. We find evi- 
dence of a first order transition line that ends in a second order point. 
By means of a rotation in parameter space we introduce thermody- 
namic magnitudes and critical exponents in close resemblance with 
simple models that show analogous critical behaviour. The measured 
data allow us to fit the critical exponents finding values in agreement 
with the mean field prediction. The location of the critical point and 
the slope of the first order line are accurately given. 



1 Introduction 

The nonperturbative formulation of four dimensional field theories is a very 
debated subject. At present, only asymptotically free theories can be rigor- 
ously constructed (See |J for a review). 

On the other hand, there is evidence @, |3| that a self-coupled (Ay? 4 ) four 
dimensional scalar field theory is trivial, that is, it describes a free theory 
after — nonperturbative — renormalization. 

The question of nontriviality of scalar fields coupled to gauge fields is, 
however, not so clear Q. In the last few years a considerable effort on the 
understanding of this problem has been carried out. 

In this work we study a four dimensional theory with a continuous sym- 
metry group: the fixed module U(l)-Higgs theory. Although it does not 
represent a limit of the SU(2)xU(l) theory (it lacks a global SU(2) sym- 
metry), we expect that many of the results obtained here may be useful for 
more complex models. 

The phase diagram of our model is represented in Fig. 1 where it can 
be seen that there are three phases: Confining, Higgs and Coulomb (strictly 
speaking only two since the first and second are analytically connected). 

We have focused our attention in the line that separates the first two 
phases. We found that it corresponds to a first order phase transition. At 
the end-point of the line we have observed a clear critical (second order) 
behaviour, where it is possible to define a continuum limit. We study the 
critical exponents at this point, as they are useful to discover the proper- 
ties of the continuum theory, in particular whether the theory is trivial or 
interacting. 

We study the model with the parameter A = oo (see (|1|)), which is equiv- 
alent to fix the modulus of the Higgs field |$| = 1. It is generally assumed 
that the fixed modulus theory belongs to the same Universality Class as the 
full (A finite) one. 

In order to find the critical exponents we compute the evolution of several 
thermodynamic quantities for different coupling and several lattice sizes. 

We have performed Monte Carlo simulations on lattices ranging from 6 4 
to 24 4 . The results presented here amount to 8 months of CPU of a custom 
64 INMOS T805 transputer machine with a performance of 100 Mflops. 

In section 2 and 3 we formulate the model and describe some of its known 
properties. Scaling relations and critical exponents near the critical point are 



also introduced. In section 4 we define the observables that will be measured 
in the simulation. The numerical method is described in section 5. Finally, 
the results are shown in section 6. We include in the appendix the description 
of the dedicated multiprocessor machine designed and constructed by our 
group. 

2 Formulation 

The action for a self-coupled scalar field with a local U(l) gauge symmetry 
(U(l)-Higgs model) in a lattice can be written as 



S = -(3J2 Ref/ r ^-«:^Re<l> r ^$ r+/ , + A^(|$ r | 2 -l) 2 + 4«:^|<l> r | 2 , 

r,fj,<u r,/i r r 

(1) 

where r is the four dimensional lattice site, the Greek indices /i, v £ {1, 2, 3, 4} 
represent the four lattice directions and $ r is the value of the complex scalar 
field at r, U T ^ is the gauge (U(l)) field at the link in the \i direction beginning 
at r, and U T ^ V is the plaquette defined by the site r and the directions /x and 
v. 

It is usually assumed that the action (|l|) for finite A belongs to the same 
universality class that the one in the A — > oo limit. Taking this limit, we can 
fix the modulus of the scalar field, and the action becomes (up to an additive 
constant) 

S = ~f3 J2 Re Ur,^ - « E Re ®rU r>IM $ r+tl (2) 

In this way both the gauge and the scalar field belong to the U(l) group. 
We will limit ourselves in this work to study this two-parameter model. 



3 Critical behaviour 

3.1 Description of the parameter space 

Let us briefly describe the different limits of the action when the coupling 
parameters /3, k take the extreme values or oo. 



3.1.1 /? = oo 

In this case the gauge fields are frozen out and the remaining model is: 

S = -/^Re5> r $ r + M (3) 

This action has a global 0(2) symmetry (or U(l)). The 0(N) model in four 
dimensions shows a continuous transition between a disorder phase, with 
explicit 0(N) symmetry at low k and a ordered one, at high k, where the 
0(N) symmetry breaks down to a O(N-l) symmetry . Due to the Goldstone 
theorem this phase has N-l massless Goldstone bosons (spin waves). This is 
a gaussian second order phase transition ||. 

We can estimate k c using the mean field approximation (MFA), the result 
is i: 

N u\ 

«c,MFA — — (4J 

Zq 

where q is the coordination number, 2d in our case (rf-dimensional square 
lattice). This is an approximation from below since MFA neglects the fluc- 
tuations and overestimates the interactions. 
An upper limit is [0]: 

* < ^ (5) 

where: 

_ /•*■ dk 1 ...dk d 1 , . 

IW ~L (2ny EMl-coBfo)) ( j 

If we take N = 2 and d = 4 we get: 

0.125 < k c < 0.311 (7) 

3.1.2 (3 = 

This situation corresponds to a spin model with annealed bond disorder [0 
(the disorder is in thermodynamic equilibrium with the Higgs field). The 
action is: 

S = ~Kj2^^ r U r ^ r+ll (8) 

r,fi 



There is not a phase transition in this case, for the following reason: it can 
be made a gauge transformation which maps the Higgs field into a constant 
because the Higgs field lives in the U(l) group. The functional integration 
over the Higgs is trivial and the U(l) sector is a one-link separable theory 
(without interaction). 

3.1.3 k = 

The remaining model is U(l) pure gauge. This model shows a transition 
between a maximally disordered Confining phase, at low (3, and a Coulomb 
phase (it has free photons and coulomb potential between static charges) at 
high (3. The transition, located at (3 ~ 1, is found to be first order ||. 

Applying a theorem by Shrock , the full theory at low k can be written 
as a pure U(l) theory with a shifted coupling, where the shift is proportional 
to k a . Thus we can extend the k = transition to a small k. 

3.1.4 k = oo 

It is always possible to fix the unitary gauge, so that the Higgs fields disap- 
pears and the action (0) becomes 

S = -(3 ]T Re U T ^ V - k J2 R e U Ttlt (9) 

Now in the k = oo limit, the second term in the action makes that 
only configurations that satisfy Re U TjlM = 1 have non vanishing probability. 
Then, S becomes trivial, and there is no transition in this case. The situation 
changes when charged (q > 1) Higgs fields are considered 0. Then the limit 
K = oo corresponds to a non trivial Z q gauge theory. 

3.1.5 Interior of the parameter space 

The phase diagram of the model was studied on the pioneering work of Frad- 
kin and Shenker [It|, and is plotted in the figure 1. Further studies can be 
found in references [p], |[ |TT|. 



On the point "C" there is coexistence of Confining, Coulomb and Higgs 
phases (triple point). 

If we consider A ^ oo the main modification is that the end point "D" 
tends, as A decreases, to the (3 — axis and finally cuts it. 
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The transition line "B-C" is a line of second order transition, the vertical 



line "A-C" is of first order [5 | and we will show that the line "D-C" is, again, 



of first order. The end point ("D") is according to our study a second order 
point. 

3.2 Critical behaviour in related systems 

Regarding the "C-D" line, the structure of our model is similar to that of a 
wide variety of systems, like, for instance, the magnetic Ising model or the 
liquid-vapour transition. 

The Ising model with nearest neighbour interaction has an action 

s = -jJ2 SiSi-hj^Si (io) 

it shows a first order phase transition line (h = 0, J > J c ) that finishes in a 
second order critical point at J c . 

The position of the end point may be determined studying the behaviour 
of the magnetization over the first order -straight — line. It is found that 

M~(J-J c f, J>J C (11) 

The exponents a and 7 are defined respectively from the critical behaviour 
of the specific heat and susceptibility: 

C~\J-J c \~ a (12) 

X~\J-J c r< (13) 

which hold the scaling relation [|l!| 

a + 2(3 + 7 = 2 (14) 

Another analogous system is the liquid-vapour transition. The lack of 
a symmetry implies that neither the straightness of the critical line nor its 
exact location (h = in the Ising case) are given. In this system the critical 
exponents (a,/3, 7) have been experimentally measured [|IIJ and the relation 
(111 checked. 



In our model we have not an explicit symmetry and so, the exact position 
and shape of the transition line, must be numerically computed. However 
we can define the critical exponents as in the previous models. 
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4 Observables 

We define the (normalized) plaquette and link energies as 



^^E^^ ( 15 ) 

EL = ^J2 Re ®rU r ,^ r+ , (16) 

where V = L A is the volume of the lattice. In terms of the above energies, 
the action can be rewritten as: 

- S = f3E P QV + kE l 4V (17) 

Both Ep and El lie in the [— 1, 1] interval. With our definition, Ep — ► 1 when 
P — > oo and El — > 1 when « — > oo. < -Ep >=< £x >= at /3 = « = 0. 
Let us write the partition function as 

Z(P,k) = j[dU}[d^]e~ s . (18) 

It is useful to introduce the parameter k' = 2k/3 to symmetrize the action. 
In this way the fluctuation matrix (or connected correlation) can be written 

as 

F„ = < E.E, >-< Ei >< E, >= ^Xt;^ < 19 > 

where 

Xl =p, X2 = K \ E 1 = E P , E 2 = E L (20) 

At a given point (Pq, k' ), the F matrix can be diagonalized. We shall call 
Amax and A min to their maximum and minimum eigenvalues respectively. 

We can perform a rotation of angle 6 in such a way that in the new 
coordinates 



cj_= /^cos^ + K'sintf (21) 

C|| = -p sin 9 + k! cos 9 (22) 

the fluctuation matrix is diagonal. 



Consequently the operators 



E ± = E 1 cos 9 + E 2 sin 9 
E\\ = -E x sin 9 + E 2 cos 9 



(23) 
(24) 



are uncorrelated. In terms of the new quantities the action can be written as 



S= (c ± E ± + a\E\\)6V 



(25) 



Assuming that the point "D" is second order, we expect divergences in 
some magnitudes. Using the Ising model analogy discussed in section (3.2) we 
can write the following formulae for the previously defined thermodynamic 
quantities: 



AE\(c\\ 



df 



dc\ 



df 



dc\ 



( c ll - c \\ 



Cll < Cm 



c^=a~ 



X{c\\) 



C(c,|) = 



d 2 f 



d& 



d 2 f 



del 



(<6V)\ X 



{<6V)\ 



-9 -7 



~ Cm 



(26) 



(27) 



(28) 



where / = -^ log Z is the intensive free energy, \ is the susceptibility and C 
is the specific heat. We have call a to the value of the c± parameter on the 
first order line that, as we shall see below, is almost independent of cii in the 
interesting region (neighbourhood of "D"). 

The critical law for AE± is analogous to that for the Ising magnetization. 
We denote as AE±_ the difference between E± in the Confining and Higgs 
phases in the c± — > a limit. 

Therefore, we expect the three critical exponents defined above to follow 
the scaling relation (|I4). 



5 The numerical method 

We have used the subgroup Z N with iV = 1024 as an approximation of 
the gauge group U(l) since, in the region of interest, the fluctuations of 
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the variables are large (typically of the order of 1 radian), and the phase 
transition associated with the discrete group is safely far away. 

The updating algorithm is an adaptive step size Metropolis, with an ac- 
ceptance rate of more than a sixty percent. 

5.1 Parallelization 

We have used lattice sizes 6 4 , 8 4 , 12 4 , 16 4 and 24 4 implemented on a transputer 
(IMS-T805) |14| machine of 64 processors with a 8 x 8 topology. For the 
smaller lattices, we have also used transputer boards of 8 and 3 processors 
with 4 x 2, 8 x 1 and 3x1 topologies. 

The parallelization strategy is straightforward: divide the lattice among 
the processors, so that each of them accounts for a smaller sublattice that 
can be updated in parallel with the other transputers. One of the problems 
is that, being a nearest neighbour interaction, the update cannot be entirely 
made inside each processor: the sites and links in the border of the sublattice 
must know some of the variables in the neighbouring transputer. 

A critical choice in the parallelization is the way in which neighbouring 
processors exchange the necessary information during the update. The sim- 
pler method could be the transmission, by means of a parallel process, of the 
variables that are required in a given step of the calculation. In addition of 
a minor problem of synchronization, transmitting single variables is poorly 
efficient regarding the link bandwidth and also, it means a considerable over- 
head due to the frequent start and end of needed parallel processes. 

For this reason we have added to the sublattice in each transputer the 
rows and columns of the links and sites needed to make the update entirely 
inside. For instance, in the 16 4 lattice on the 64-processor machine with 
a 8 x 8 topology, each transputer holds a sublattice of 4 x 4 x 16 2 lattice 
and not a 2 x 2 x 16 2 lattice. The update is performed almost synchronized 
in all processors. Once it is finished, the borders are transmitted to the 
neighbouring processors. To guarantee the independence of the variables 
currently updating, we perform a checkered update. 

We have obtained in this way a parallel efficiency close to a 95%. 
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5.2 Spectral Density Method 

The precise location of the critical values of the parameters may be a difficult 
task because Monte Carlo methods provide information about the thermo- 
dynamic quantities only at particular values of the couplings. The approach 
that we use here to locate them is based on histograms and is known as the 
Spectral Density Method flT5|| . 

The generalization of the method to a two-dimensional parameter space 
is straightforward. Considering that, we perform a Monte Carlo simulation 
at a particular point of the parameter space (/?,«), and we compute the 
histogram H(Ep, El) as an approximation to the density of states. The 
probability of finding a configuration of plaquette energy Ep and link energy 
El at a different point (/3i , K\ ) can be calculated as 



H ( E E ) (Pi-P)6VE P +( K1 - K )4VE L 

1 P ' L) fdEpdE L H{Ep,E L )e^-W VE 'r + ^-^ VE, L' { ] 

Let us discuss the range of applicability of the Spectral Density Method. 
Let ap, ol be the widths of the measured histogram in the Ep, E L directions 
respectively. It may be easily seen from the previous equation that the ranges 
are A/3 ~ 1/(6W P ), A/s ~ l/(4Va L ). 

Although near the transition line op (or analogously ol) is large, the 
application of (|29D at fixed k is very useful to find the /3-value where the 
fluctuation of the energy has a maximum — apparent critical point — and, 
eventually, to adjust the parameters for a new simulation. 

We can also move simultaneously in both directions. The minimum eigen- 
value of the fluctuation matrix corresponds to an eigenvector parallel to the 
transition line. This means that the range of applicability of the Spectral 
Density Method in the c\\ direction is large. 

We can go one step further using data from simulations at different points 
of the parameter space, using a two-dimensional generalization of the mul- 



tihistogram method proposed in [ 16] , which gathers all the information for a 
given lattice size. This method has been used for many observables and the 
results plotted as a smooth line in the figures, with the points corresponding 
to single simulations. However, in order to obtain a safer estimation of the 
errors, we only use the single points to do the fits to a critical power law. 
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Usually, the Spectral Density Method is very useful to find the value of 
the coupling where some observable has a maximum, as well as to obtain an 
accurate value for this maximum. In this work we have used it extensively 
to locate the transition line "C-D" at a fixed en, looking for the maximum in 
the fluctuation of E±. However, in order to find the end point "D" we have 
to move in the cm direction, in doing so we do not find a maximum for any 
simple observable, and the usefulness of the method is a great extent lost. 

5.3 Measurements 

At every simulation point we store the plaquette and link energies to con- 
struct the histograms. We compute the two-dimensional histogram in the 
Ep-El plane. As an example, we show in figure 2 some contour plots of the 
energies histogram. 

In principle we can move in the (j3, k) plane in whatever direction using 
the spectral density method. In figure 3 we plot the mean value of E± as a 
function of the couplings obtained with the multihistogram method. 

Nevertheless, most of the results presented in this work, have been ob- 
tained studying one dimensional histograms. To this end, we rotate the 
parameter space in order to discretize the energies E±. To simplify the com- 
putations we have chosen a fixed rotation angle for all lattice sizes and all 
the parameter space points. Notice that the error in the determination of 
the angle could only mean second order corrections for most of the quanti- 
ties of interest. The only exception is the minimum eigenvalue A m j n that is 
computed diagonalizing the fluctuation matrix in every simulation. 

We usually discretize the energy interval into one hundred subintervals. 
From the (one dimensional) histogram (figure 4, upper side), we compute 
E± and dE±/dc± in the neighbourhood of the simulation point, determining 
the best approximation to the critical point by looking at the maximum of 
the derivative. In figure 5 we show an example of this method. The validity 
range is estimated from the fluctuation in the energy. 

A very important quantity to be computed is the latent heat, that is, 
the difference between the values of E± on both sides of the first order line. 
In a finite lattice this limit is not well defined; we take as its definition 
the distance between the two minima of the effective potential at the first 
order line (that is, the distance between the two maxima of the histogram). 
Unfortunately, this procedure needs a precise estimation of the local maxima 
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of a noisy function. To reduce the statistical error, we smooth the histogram 
in the region near each maximum with a cubic spline, measuring the distance 
between the smoothed functions (see figure 4, lower side). 

The computation of statistical errors has been carried out using the jack- 
knife method. We perform a primary determination of the correlation in 
Monte Carlo time and construct statistically independent bins. The number 
of iterations performed in the larger lattices for k G [0.52, 0.54] has been (in 
thousands of Monte Carlo Sweeps) 
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(30) 



so that, in the neighbourhood of the critical point, we can use bins of around 
a hundred thousand of MC sweeps. 

6 Results 

In this section we report our results classified according to every observable 
computed. In all cases we have used the coordinate transformation to the 
(c_[_, C||) plane. 

For the sake of simplicity regarding the figures we fix from the beginning 



9 = 0.96 



(31) 



which is within a 1% our estimation for the angle of the fluctuation matrix 
defined above. We will discuss in every case the effects of the particular 
selection of this quantity. 

We have always run in points of the parameter space over the first order 
line for c\\ < cjj and on the prolongation (dotted line in figure 1) for cy > cjj. 

Since there is not an explicit symmetry, as in the liquid vapour transition, 
the straightness of the transition line is not implied, and then the angle of 
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the eigenvector of the fluctuation matrix does not have to take the same 
value as the one of the first order line. However we have found numerically 
that, over the line, c± is practically constant; it changes less than 0.01% for 
K E [0.52,0.54] ( C || e [-0.501,-0.478]). 

6.1 Latent Heat 

As we saw in ( p^ ) the critical exponent (3 is related with the behaviour of 
AE± over the critical line. 

At the first order line there is a discontinuity in E± and the gap between 
its values on both sides, is, in the thermodynamic limit, the latent heat. 
In a finite lattice it is not easy to measure this gap since the discontinuity 
is rounded. The method used in this work is to compute the gap from 
the histogram in the quantity E±. More precisely we compute the distance 
between the two maxima of the histogram. If the height of both is not 
the same we shift (with the Spectral Density Method) the histogram the 
necessary amount in c±. 

Since E\\ is continuous, it is equivalent, regarding the critical behaviour, 
to make the analysis in terms of E± or almost any linear combination of E± 
and E\\, with the only requirement of a nonvanishing coefficient in E±. In [|17| 



we used the latent heat for the plaquette energy. We expect a — slightly- 
better measurement when choosing the optimum combination of Ep and El 
which is E±. 

For the same reason, a small error in the determination of the rotation 
angle 9 is unimportant. We remark that this error affects only quadratically 
so that it is completely negligible. 

In figure 6 we show our measurements for AE± for several couplings and 
lattice sizes as a function of cy . 

Although the critical behaviour is very clear, a precise determination of 
the critical exponent (/3) is very difficult. 



Equation ( 26|) is only followed strictly in the thermodynamic limit. In 
a finite lattice, we should find deviations from the functional form E± = 
A{c\\ — cjj)^, however, it is expected that the main deviation can be considered 
as a finite size dependence on the parameters {A, c|,/3}. In particular we 
expect the strongest dependence on the parameter cjj (see section below on 
Finite Size Scaling). The procedure we use to compute the parameters is the 
following: we fit the whole data to the function fl26|) with j3 independent of 
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L, allowing a finite size dependence for cjj and A. Successively we discard 
the data from the smaller lattices to check the asymptotic behaviour. Our 
results are 

£ = 6,8,12,16 = 0.53(7) 

L = 8,12,16 = 0.53(7) , , 

L = 12,16 = 0.55(9) l > 

L=16 = 0.50(11) 

Although the statistical errors do not allow us to observe a monotonous 
evolution to the thermodynamic limit, the previous results give a strong 
evidence in favour of the classical value = 1/2. 



In reference [[17] we compute the exponent using only the data from 
the plaquette energy. Although the simulations are essentially the same, 
some minor variations are found since the observables are not completely 



correlated. In g7| we found = 0.54(6) for L = 12, 16 and = 0.47(9) for 
L = 16. 

In figure 7 we plot the latent heat squared as a function of the parameter 
Cm for the points near the critical one. The linear behaviour seems to be in 
agreement with the data. 

6.2 Maximum eigenvalue 

An alternative way to compute the exponent is measuring the fluctuation 
of the energy. The fluctuation of E± is what we previously called A max . In 
the limit of a histogram with infinitely narrow peaks, A max = (Ai5j_) 2 /4. 
We point out that measuring A max instead of AE±_ is similar to measure the 
square of the magnetization (in a magnetic system) to avoid the cancellation 
of the magnetization due to tunneling effects. 

In practice it is not necessary to diagonalize the fluctuation matrix at 
each point since the variations of the angle 9 only affect quadratically, and 
mainly as a global multiplicative constant. 

In figure 8 we show the evolution of A max as a function of c\\ for several 
lattice sizes. We see a window in the larger lattices where the behaviour is 
almost linear according to A max = A(c\\ — c Tj) 2 with near 1/2. 

We remark that the statistical error in the measure of A ma x is much smaller 
than the one for the latent heat (compare figures 7 and 8). Unfortunately, 
the deviation from the square of the latent heat due to finite size effects is 
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large, so it is not easy to obtain a precise estimation of the exponent (3 and 
of its error. 

6.3 Susceptibility 

The susceptibility x = dE±/dc± is related to the maximum eigenvalue of 
the fluctuation matrix in the absence of phase coexistence and there its mea- 
surement is straightforward: x — 6V^m&x- When c\\ < cjj we would have to 
measure x a t c± near the transition line and, after that, take the limit from 
one side of the transition. Alternatively, and this is the method we use, we 
can divide the histogram in E± in two halves and measure the fluctuation on 
each half. In the thermodynamic limit the results are equivalent. However, 
for finite lattices near the critical point the overlap between both peaks is 
big and the measure of x cannot be very precise. 

In figure 9 we show our results for 6VX max which can be called suscepti- 
bility only for cy > cjj. 

The finite size effects do not allow us to use directly the points in the 
simulated lattices for fitting to the critical power law (see figure 9). To 
compute 7 we first take the thermodynamic limit and then fit the asymptotic 
values. We must point out that this indirect process reduces the objectivity 
in the determination of the error. Our results are 

7 = 1.13(17) (33) 

according with the Mean Field prediction 7 = 1. 

In figure 10 we plot the results for x i n the larger lattices, including also 
the values for c\\ < cjj obtained using the division method discussed above. 
We can obtain a more precise value computing the dispersion excluding the 
region of the histogram between peaks: we compute the dispersion in the 
first phase integrating over the left part of the first peak, and in the second, 
integrating over the right part of the second peak. The accumulated error 
is difficult to compute, but we can obtain a crude estimation comparing the 
results from both peaks (see figure 10). 

We stress that with the hypothesis of the same exponents on both sides 
of the transition (7 = 7') we obtain again (after excluding the points closer 
to the critical point which show great finite size effects) a value for 7 near 1. 
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Finally we remark that the divergence of the susceptibility when ap- 
proaching the thermodynamic limit is a clear evidence of a second order 
behaviour at the end point "D" . 

6.4 Specific Heat 

The small eigenvalue of the fluctuation matrix is related with the energy E» : 

1 dE\\ 

Amin — 777"^ (34 J 

6V oc\\ 

so that, in our analogy with the Ising model, we can call it Specific Heat. 
Since dE\\/dc\\ ~ A|c|| — c fjl~ a > V^mm should present a divergence at the 
critical point if a > 0. If a — the divergence should be logarithmic and so, 
hard to observe. 

To determine the small eigenvalue of a matrix with a much larger one is 
not an easy task. However we have obtained quite precise determinations. 
Notice that the minimum eigenvalue is just the width of the two-dimensional 
histogram (see figure 2) which can be clearly distinguished from its length. 

In figure 11 we plot the minimum eigenvalue (times V) as a function of 
cm. The absence of a divergence at the critical point practically excludes the 
possibility of a positive value for a. 

In a strict sense, when there is phase coexistence we would have to com- 
pute the minimal eigenvalue for both phases. For large lattices the Ep — El 
histogram becomes very narrow, and the difference in the angle at each phase 
makes the minimum eigenvalue for the whole histogram to grow (see figure 
11). Nevertheless, this does not change our conclusions about the behaviour 
at the critical point. 

6.5 Finite Size Scaling 

The critical exponent v has been computed with a Finite Size Scaling Anal- 
ysis. We have study the shift in the apparent critical point as a function of 
the lattice size. We expect that the shift will follow the law 

Acfi(L) ~ L~« (35) 
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In figure 12 we plot the values (squared) obtained fitting the latent heat 
data for each lattice size fixing the value of j3. A least squares three parameter 
fit gives 

v = 0.52(4) (36) 

The value (^) is almost insensitive to the value of /3 used in the extrapolation. 
For values of (3 in the range [0.4, 0.6] the resulting v varies by less than 0.01. 
In reference JL7| using just the plaquette energy data, we obtained the 
value v = 0.47(4). 



7 Conclusions 

We have found a second order point where critical exponents a, (3, 7 may be 
defined in close analogy with ferromagnetic spin systems. The location of 
this point has been accurately measured: 

c\ =0.77391(2) 

cjj =-0.494(1) [6i) 

it corresponds to (3 C = 0.8485(8), n c = 0.5260(9). Although we have fixed in 



(pTD the rotation angle 6 that defines c± and cm, our best estimation of the 
angle that diagonalizes the energy fluctuation matrix at the critical point is 
e c = 0.963(3). 

The thermodynamic limit extrapolations of the latent heat data corre- 
sponding to simulations performed below cjj clearly give nonvanishing values, 
that show the first order nature of the critical line between the Higgs and 
the Confining phase. 

Fitting the thermodynamic limit extrapolations of the inverse suscepti- 
bility we have seen that it approaches zero at cjj, giving us a strong evidence 
of its second order behaviour. This result has also been confirmed with ex- 
trapolations from the region below cjj. 

Our estimations of the critical exponents are compatible with the mean 
field results a = 0,(3 = 1/2, 7 = 1, v = 1/2. 
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8 Appendix 

The 64-processors machine that we have call Reconfigurable Transputer Net- 
work (RTN), has been entirely designed and build inside our group. It pro- 
vides us with a power of 100 (sustained) Mflops at a low cost. In the actual 
configuration it holds eight boards with eight T800 each, plus a controller 
board with one transputer. RTN interfaces with a host computer, a PC in 
the actual configuration, via another board (root board). The eight trans- 
puter boards have also a cross-link C004 (programmable switch) that allows 
to interconnect them. Each transputer has 1 Mb of memory (may have up to 
4Mb). They are connected inside the board forming a ring. One of the two 
remaining links of each transputer is connected to the C004 and the other 
free link is connected to the C004 of the next board. In this way we may 
have a torus (8 x 8) topology for the 64 transputers. Nevertheless it is pos- 
sible through the C004s to attain different topologies. We could divide, for 
instance, RTN in eight identical machines dedicated to different problems. 

The purpose of the controller board, which has a transputer and a cross- 
link , is to boot the code to each transputer of RTN, close the torus, wait for 
the end of the calculation, open the torus, read the results and repeat the 
process. The way to know when the calculation is done and when to open 
the torus is through the event pin of one of the transputer in each board. 

The root board, as we said, interfaces with the host computer and is 
responsible for sending the code and input to the controller board, it also 
reads the output from the controller and writes it on the host devices (disk, 
terminal, etc.). 
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Figure captions. 

1. Scheme of the phase diagram of the fixed module U(l)-Higgs model. 

2. Contour plots of two-dimensional (Ep-E^) energy histograms for L = 8 
(upper part) and L — 12 (lower part). The left side ones correspond 
to the point ((3, k) = (0.854, 0.52) and those on the right to (/3, k) = 
(0.835, 0.54). In all cases the most external contour is plotted at a value 
of a 10% of the maximum. 

3. Ej_ near the critical point in a L = 8 lattice using a multihistogram 
method. 

4. Example of one-dimensional histogram (above). Example of the mea- 
sure of the latent heat with a cubic spline fit (below). The data corre- 
spond to 500000 sweeps on a L = 16 lattice for (f3, «) = (0.8495, 0.525). 

5. E± (above) and dEj_/dc± (below) using the Spectral Density Method 
obtained from a simulation with L — 12, /3 — 0.84945, k = 0.525, 
500000 MC sweeps. The filled circle is plotted at the cm of the simula- 
tion. 

6. Jump of E± across the first order line. The continuous lines are ob- 
tained from a power law fit. 

7. Square of the Latent Heat as function of cm. The straight lines are 
minimum squares fits. Symbols as in figure 6. 

8. Maximum eigenvalue of the energy fluctuation matrix against c\\. The 
continuous lines are obtained with a multihistogram method. 

9. Inverse of the susceptibility (when there is no phase coexistence) as a 
function of cn . The continuous lines are obtained with a multihistogram 
method. 

10. Inverse of the susceptibility on both sides of cjj . Only represented for 
L = 16 (circles) and L = 24 (triangles). 

11. Minimum eigenvalue of the energy fluctuation matrix against cm. 

12. cjj(L) as a function of 1/L 
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